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. A significant fraction of the energy emitted in the early Universe came from 
CN ■ 

^ [ very luminous galaxies that are largely hidden at optical wavelengths (because 

2 ■ of interstellar dust grains); this energy now forms part of the cosmic back- 

^ ■ ground radiation at wavelengths near Imm.^ These submillimetre (submm) 

Q ! galaxies were resolved from the background in 1997,^ but have been difficult 
to identify and study due to the poor spatial resolution of submm instruments."^ 



This has impeded the determination of their distances (or redshifts, z), a cru- 



^ , cial element in understanding their nature and evolution. Here we report 
I spectroscopic redshifts for ten representative submm galaxies that we identi- 
fied reliably using high resolution radio observations. ^'^'^ The median redshift 
for our sample is 2.4, with a quartile range oi z = 1.9—2.8. The submm popu- 
lation therefore coexists with the peak activity of quasars, which are thought 
to be massive black holes in the process of accreting matter, suggesting a close 
relationship between the growth of massive black holes and luminous dusty 
galaxies.^ The space density of submm galaxies at 2 > 2 is about 1000 times 
greater than that of similarly luminous galaxies in the present-day Universe, 
so they represent an important component of star formation at high redshifts. 
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The first real insight in the origin of the far-IR/submm background came with the 
commissioning of the SCUBA submm camera on the JCMT. Almost 100% of the submm 
background has been resolved with the deepest SCUBA maps, which exploit a sensitivity 
boost from gravitational lenses. ^''^'^ However, the near-IR/optical faintness of the submm 
galaxy population conspires with the modest positional precision available from the large 
(15") SCUBA beam and the large surface density of unrelated optically-faint galaxies to 
render positional coincidence alone inadequate to identify counterparts.^'^ This problem 
is compounded by the small field of view of SCUBA: compiling large samples of submm 
galaxies is slow, with roughly one detected every night. Because of these difficulties, 
robust spectroscopic redshifts have been published for only a handful of submm sources, 
all atypically bright at optical wavelengths. 

Here, we have overcome these problems by taking advantage of 1.4-GHz radio data 
from the Very Large Array (VLA) radio telescope (Hubble Deep Field^^, SSA13 field^^, 
ElaisN2 and Lockman fields^). These deep images have fine (1-2") spatial resolution and a 
large (30' FWHM) field of view. The radio data was exploited to foUowup SCUBA submm 
data in the Hubble Deep Field (12 hr RA) and SSA13 (13 hr RA) fields^5'^'^^ and the 
Lockman (10 hr RA) and ElaisN2 (16 hr RA) fields^^. The radio data are sensitive to 
the synchrotron radio emission from cosmic-ray electrons accelerated in the supernovae 
explosions of the same high-mass stars that heat the dust and generate the far-IR/submm 
emission. Hence a deep radio image should be an efficient route to pinpoint the positions 
of many submm galaxies. This is confirmed in fields with both submm and radio observa- 
tions: radio counterparts brighter than 30 fiJy can be found for 60-70% of submm galaxies 
brighter than 5mJy at 850 /xm (which have a surface density of 450deg~^,^ and contribute 
about 20% of the submm background). ^'-^^ Moreover, the accurate radio positions mean 
that scuba's efficient 'photometry' mode can be used to search for submm galaxies, 
raising the detection rate to around ten per night. 

Although submm galaxies have faint optical continuum emission,^ the radio-detected 
samples are sufficiently large to allow efficient, multi-object optical spectroscopy using 10-m 
telescopes. Slit positions can be assigned accurately, as the radio/submm emission can be 
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located on optical images to within 0.3-0. 8" 7 Using this approach we targetted 34 radio- 
detected submm galaxies mostly brighter than 5 mJy at 850 /xm with optical magnitudes 
in the range / = 22.2-26.4 using the ESI and LRIS spectrographs on the Keck telescopes 
(Fig. 1). We measured redshifts for ten submm galaxies, representative of the blank-field 
population. We have doubled the number of optically bright (I < 23.5) submm galaxies 
with accurate redshifts. Our six redshifts for / > 23.5 submm galaxies represent the first 
redshifts for this class of submm galaxy. The redshifts span the range z = 0.8-3.7, with 
a median of 2.4 and an interquartile range of 1.9-2.8 (Fig. 2). These redshifts allow dust 
temperatures Ta and bolometric luminosities to be determined (see Table 1), assuming 
that the tight correlation between far-IR and radio emission observed at low redshifts ""^^ 
remains valid. Typical values are of order 35 K and several 10^^ respectively. Note 
that relativistic electrons accelerated in shocks close to an active galactic nucleus (AGN) 
should boost the radio flux density above that expected from the standard correlation. 
If the radio emission is boosted by an AGN then both temperature and luminosity are 
overestimated. There are no obvious cases in the present sample in this category. 

The redshifts impose robust limits to the space densities of submm galaxies in the 
redshift ranges 0.5-1.2, 1.8-2.8 and 2.8-4 (the gap at z = 1.2-1.8 is the spectroscopic 
desert in which no strong restframe UV lines are redshifted into the optical window): we 
calculate p = (3.3 ± 2.3) x 10"^ Mpc'^, (6.5 ± 2.5) x 10"^ Mpc-^ and (2.4 ± 1.2) x 10"^ 
Mpc~"^, respectively, assuming a flat A = 0.7 Universe with Hq = 65 km Mpc~^. This 
compares to a space density of ~ 10~^ Mpc""^ for low-redshift galaxies with comparable 
bolometric luminosities greater than 5 x 10^^ Lq, comparable to the luminosities of our 
spectroscopically targeted galaxies at any redshift greater than one. The volume density 
of very luminous galaxies thus increases almost a thousandfold from 2; ~ to 2. 

The relationship between different classes of high-redshift galaxy is critical for our 
understanding of galaxy evolution. Our spectroscopy shows that the submm galaxies are 
coeval with the important populations of star-forming galaxies and quasars detected at z = 
2-3 using optical and X-ray techniques. ^'^^'^^'^^ We can use their relative space densities to 
compare and relate these populations. In a 1000 Mpc^ box at 2; ~ 2.5, there are 10 Lyman- 
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break galaxies with Rab < 25.5 (~ 0.2xL^5qq) and 1 submm galaxy from our sample (with 
850/xm flux >5mJy), but the single submm galaxy produces a comparable luminosity to 
the 10 Lyman-break galaxies. (Here -Rab is the flux normalized magnitude near 6900A; 
-^1500 is the characteristic luminosity at rcstframc 1500 A, denoting a change in slope in the 
luminosity function.) Equally, the volume density of submm galaxies is similar to that of 
narrow-line Type-II AGN selected at these redshifts, from both X-ray and optical surveys, 
and an order of magnitude greater than that of bright optical quasars with Mb > —25.^ 
However, the links between the different classes of sources all depend critically on whether 
the activity is intermittent or continuous. 

To quantify the whole submm population using this survey we must account for several 
selection effects. In particular, requiring a radio identiflcation prior to spectroscopy limits 
the maximum redshift. Studies of galaxies at low and moderate redshifts suggest that 
changes in the dust properties, especially the dust temperature Tj, modify the relative 
strength of the emission in the submm and radio wavebands. The detectability of distant 
submm galaxies (at a flxed 850-//m flux) in the radio waveband thus depends on the dust 
temperature: warmer galaxies are detectable at ^ > 3, while the radio emission from cooler 
sources falls below the flux limit at lower redshifts. About 30% of submm galaxies that 
are brighter than 5 mjy at 850/im are not detected in the radio: these galaxies would be 
lost from our catalogue at ~ 2.5 for the coldest plausible ~ 25 K, and at 2; > 3 for 
the typical ~ 35 K that we flnd. Hence, we can only place a lower limit on the space 
density in the highest redshift interval. The density derived ai z = 2 should be unaffected 
by the need for the galaxies to be detected at radio wavelengths. 

We must also consider the completeness of our spectroscopic observations. While our 
redshift completeness is relatively high (~30%), we failed to flnd convincing redshifts for 
24 galaxies. However, nine of of these have strong single emission line detections which 
could plausibly be identifled as Lyman-a. Some of the remainder without emission lines 
could lie in the spectroscopic desert, making redshifts difficult to identify. The absence of 
sources at ^ ~ 1.5 in Table 1 suggests that this effect is probably responsible for some of 
the incompleteness. Our targets generally have faint UV continua, and so there is a bias 
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towards identifying redshifts for galaxies with strong emission lines. If the missing galaxies 
lie outside of the desert, but lack strong emission lines, then their radio detection suggests 
they probably lie in the same z — 1-4 range spanned by the spectroscopically identified 
sources. 

On the basis of bolometric luminosities estimated from our radio and submm mea- 
surements and redshifts (Table 1), assuming the far-IR-radio correlation, we calculate the 
following luminosity densities from our survey: Pl(0.5 < ^ < 1.2) = (5±5) x 10^ Lq Mpc~^, 
Pl(1.8 <z< 2.8) = (7i|)xlO^L0 Mpc-^ andpL(2.8 <z< 4.0) = \.q Mpc-^. 

To translate these values into star-formation densities we must subtract any contribution 
to the luminosities of these galaxies from the heating of dust by an active galactic nucleus 
(AGN). The proportion of the submm population detected in deep Chandra and XMM- 
Newton X-ray surveys, which are sensitive to even dust-obscured AGN,'^'^° at flux densities 
significantly greater than that expected from star formation alone is at most 30%. Hence, 
we assume conservatively that the far-IR luminosity is derived from star formation in the 
remaining 70%. Luminosity-density measurements can be translated into a star-formation 
density using the standard calibration of 1.6-2.2x10^ (M© yr~^)~^.^^ 

The bright radio-detected submm galaxies presented here represent just 20% of the 
850- //m background, yet their estimated star-formation densities at redshifts of 2 and 3 
are comparable with those inferred from optical observations. ^-"^ Both of these estimates 
require corrections, in the case of the optical surveys an increase of a factor of 5 is needed 
to account for dust obscuration,^^ while our submm estimate needs to be corrected by a 
similar factor to account for spectroscopic incompleteness and sources below our submm 
fiux limit. We find that the submm and Lyman- break populations encompass all of the star 
formation activity needed to reproduce the star formation history inferred from models of 
the background radiation and submm counts.^ Hence, our substantially complete redshift 
survey provides the strongest evidence for a rapid increase in the volume density of dusty, 
very luminous galaxies in the distant Universe. Their density peaks at a redshift 2; ~ 2.4 
at a value about 1000 times greater than that at the present epoch. Thus we conclude that 
the balance between the star-formation density in the most luminous dust-obscured versus 
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unobscured galaxies has altered radically over the last 80% of the history of the Universe. 
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Table 1. Positions of submm galaxies with redshifts from our survey, along with their measured red- 
shifts, submm fluxes and optical magnitudes, infered dust temperatures and bolometric luminosities 
-Lbol- and -Lbol were calculated, based on the redshift, submm and radio flux densities, assuming 
that the far-IR-radio correlation applies. Sources are in order of decreasing flux density, as in Figure 1. 
(SMM163704. 3+410530 was previously published in .) Submm fluxes for sources in the Lockman (10 hr 

1 7 

RA) and ElaisN2 (16 hr RA) fields are from while fluxes for other sources were measured ourselves. Six 

out of 10 galaxies have /-band magnitudes fainter than 23.5, consistent with the optical flux distribution 

16 

of blank-field submm galaxies presented elsewhere, and confirming that our sample is representative 
of the population. Galaxies with hotter inferred dust temperatures could have their radio flux densities 
enhanced by contribution from an active galactic nucleus. Understanding the relationship between submm 
galaxies and other high-redshift populations will be helped by comparing the strength of the clustering, 
both within and between the different classes. This will require substantially larger samples than are 
currently available. 

RA dec redshift SgsO/im (mJy) /-mag (K) Lbol (IO^^Lq) 
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Figure Legends. 

Figure 1. The ten spectra from Table 1 with identified redshifts (smoothed to the instrumental reso- 
lution) from the sample of 34 submm galaxies observed, spanning the / = 22.2-26.4 range. (Note that 

27 29 
SMM163704. 3+410530 was previously published in .) Pilot observations, using Keek's ESI spectro- 

1 7 

graph on 2001 July 16 in the ELAIS-N2 field yielded three spectra with strong emission lines. Multi- 

28 

object Keck observations of 31 submm galaxies using LRIS were made on 2002 March 18—19 in light 

cirrus conditions (with good O.S" seeing) in the HDF &: SSA13 fields"^'^^, and the Lockman Hole and 
1 7 

ELAIS-N2 fields . All the spectra probe the observed wavelength range from 0.3-1 /Xm. Exposures times 
were 1.5-4.5 hr, split into 30-min integrations. Data reduction followed standard techniques using custom 
IRAF scripts. One-dimensional spectra were extracted and compared with template spectra and emission- 
line catalogues to identify redshifts. All identifications are based on multiple lines, most prominently the 
LyO; line which varies tremendously in both flux (ranging from 1 to 60 //Jy) and rest-frame equivalent 
width (3 to >100A). Weaker stellar /interstellar /AGN features and/or continuum breaks were detected in 
the spectra, strengthening the redshift identifications. Two galaxies at the same redshifts were identified in 
the vicinity of the 10 robust submm sources, as shown for the companions SMMJ105155.7+572312-a (one 
of the few submm galaxies showing Lyman-o; in absorption) and SMMJ105155.7+572312-b. The highest 
S/N spectrum (SMMJ163650. 0+405733 — ) includes features indicative of both strong starburst activity 
(P-Cygni wind absorption profiles, stellar/interstellar absorption lines) and a weak AGN, a mix typical of 
hybrid starburst/Seyfert-2 galaxies. None of our spectra have broad lines like Type-I AGN. Narrow-line 
Type-II AGN with enhanced NV and/or CIV emission are consistent with half of our spectra. Despite 
this evidence for AGN in these galaxies, they appear to be relatively weak in energetic terms as compared 
to those identified in X-ray or optical surveys. 

Figure 2. The redshift histogram of our submm galaxy sample. We describe the selection effects 
using an evolving model of the local far-IR luminosity function, in which the dusty galaxies are represented 

by a range of template spectral energy distributions. The model has been tuned to fit the statistical 

2(3 

properties of the submm-radio galaxy population (surface density on the sky and submm— radio colours). 
As a range of spectral energy distributions are represented at each galaxy luminosity, a broader redshift 
distribution results than if the luminosity was tied one-to-one with the dust temperature. We plot the 
predicted redshift distributions for submm galaxies with flux densities SgdOf^m > 5 mJy (solid line) and 
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radio sources with 30 //Jy qHz < 500 fiJy (dashed line). We expect to miss sources lying between 

the submm and radio model curves due to our requirement of a radio detection to pinpoint the submm 
source. The apparent deficit between model and data at 2 ~ 1.5 may be indicative of the difficulty in 
obtaining spectra in the spectroscopic desert that spans the range 1.2 < Z < 1.8. Note that the models 
are in good overall agreement with the observed redshift distribution. 
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SMMJ163658. 8 + 405733 z=1.189 
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